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Abstract We investigate effects of a global magnetic field on the dynamics of an ensemble of 
clumps within a magnetized advection-dominated accretion flow by neglecting interactions 
between the clumps and then solving the collisionless Boltzman equation. In the strong¬ 
coupling limit, in which the averaged radial and the rotational velocities of the clumps follow 
the ADAF dynamics, the averaged radial velocity square of the clumps is calculated analyt¬ 
ically for different magnetic field configurations. The value of the averaged radial velocity 
square of the clumps increases with increasing the strength of the radial or vertical compo¬ 
nents of the magnetic held. But a purely toroidal magnetic held geometry leads to a reduction 
of the value of the averaged radial velocity square of the clumps at the inner parts with in¬ 
creasing the strength of this component. Moreover, dynamics of the clumps strongly depends 
on the amount of the advected energy so that the value of the averaged radial velocity square 
of the clumps increases in the presence of a global magnetic field as the flow becomes more 
advective. 

1 INTRODUCTION 

Accretion processes have been extensively studied during last decades and several types of the accretion 
models have been proposed to explain certain observational features of some of the astrophysical objects. 
Most of the models assume that the accretion process occurs as one-component gaseous fluid. But there 
are strong observational and theoretical arguments which imply at least some of accreting systems are 
clumpy so that they consist of cool clumps embedded in a much hotter and more tenuous gaseous fluid. 
For example, observational evidences show that broad-line region of active galactic nuclei (AGN) has a 
clumpy structure (Rees 1987; Krolik & Begelman 1988; Nenkova et al. 2002; Risaliti et al. 201 1; Torricelli- 
Ciamponi et al. 2014 ). The broad emission lines in the spectrum of AGNs are attributed to an assembly of 
clouds which are moving through a hot intercloud medium. Clouds’ basic properties are estimated according 
to the photoionisation models. These models predict that clouds’ typical size is lO^^^^cm and their number 
density is 10^°^^cm“^ (e.g., Rees 1987; Krause et al. 2012). Orbital motion of BLR clouds is a rich source 
of information for estimating the mass of the central black hole (e.g.,Netzer & Marziani 2010). One can 
neglect collisions between the clumps and investigate the orbit of an individual clump in the presence of 
the central gravitational force and possible radiation held like a two-body classical problem (e.g., Netzer 
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there are theoretical concerns about the stability of the clumps, it is generally believed that magnetic fields 
provide a confinement mechanism (e.g., Rees 1987). 

Another approach to study dynamics of the clumps embedded in a hot medium is based on analyzing the 
collisionless Boltzmann equation as has been done by Wang et al. 2012 (hereafter WCL). They described 
the gaseous ambient medium using the classical similarity solutions of Advection-Dominated Accretion 
Flows (ADAFs) presented by Narayan & Yi (1994) for the non-magnetized systems. Although collisions 
between the clumps have been neglected for simplicity, their interactions with the surrounding gaseous 
medium was included through a drag force as a function of the relative velocity of the clumps and the gas. 
In the strong-coupling limit, it was shown that the root of averaged radial velocity square of the clumps 
is much larger than radial velocity of the gas flow. The analysis has been extended to the magnetized case 
by Khajenabi et al. (2014) where the authors considered a purely toroidal magnetic field geometry for the 
gaseous component. They found that when magnetic pressure is less than the gas pressure, the averaged 
radial velocity of clumps decreases at the inner regions of the system whereas it increases at the outer parts, 
though this enhancement is not very significant unless the system becomes magnetically dominated. 

In this work, we extend the analysis by Khajenabi et al. (2014) to include all three components of the 
magnetic field in the gaseous component. Since properties of the gas flow is significantly modified in the 
presence of a global magnetic held (Zhang & Dai 2008), the drag force varies depending on the strength 
of the magnetic held and the assumed magnetic geometry which will eventually lead to a considerable 
modification in the velocity dispersion of clumps. In the next section, we present our basic assumptions and 
the equations. In section 3, a parameter study for the dynamics of the clumps are presented. We conclude 
with a summary of the results in section 4. 

2 GENERAL FORMULATION 

Our analysis for describing an ensemble of clumps is based on WCL approach which implements the 
collisionless Boltzman equation in the cylindrical coordinates (r, <f>, z) including the components of the 
drag force. If we assume the distribution function of clumps is represented by T, then Boltzman equation 
is written as 
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where f = Vr, 4> = v^/r and z = v^- The central object with mass M is at the origin and the gravitational 
potential becomes fl) = GM/{T + The components of the drag force are = fr{vr — V).)^ and 

F<j, = where fr and are constants of order unity. Here, the radial and the rotational velocities 

of the ADAF where clumps are moving within it are denoted by Vr and V^. Thus, dynamical properties of 
the background medium affect dynamics of the clumps through these components of the velocity. Thus, 
Boltzman equation becomes (WCL) 
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It is very unlikely to solve this equation analytically in a general case unless we apply further simplifying 
assumptions. We assume that the clumps are strongly coupled with the background gaseous medium which 
implies the mean radial and the rotational velocities of the clumps are equal to the radial and the rotational 
velocity of ADAF. Under these simplifying assumptions, it is possible to obtain the root mean radial velocity 
square of the clumps analytically (WCL): 

_l-e^g-rrwj 


= c^{^[a^c?rrA3(rr,r) + (l-c^)A|(r^,r)-^^^A7(r^,r)] 


where F^ = frRsch and F^ = ftpRsch are the coefficients of the drag force. Function Ag is introduced 
by WCL as A^ = x~'^exp{—Tjix)dx. The outer boundary condition is at r = Tout, so that (u^) = 
Moreover, properties of the ADAF are described using a set of radially self-similar solutions where 
Cl and C 2 are coefficients of the radial and the rotational velocities of the ADAF. In WCL, the standard 
nonmagnetized ADAF solutions (Narayan & Yi 1994) have been used for their analysis. Then, Khajenabi 
et al. (2014) extended the analysis by including purely toroidal component of the magnetic field using 
similarity solutions of Akizuki & Fukue (2006). But we extend previous studies by considering all three 
components of the magnetic field using similarity solutions of Zhang & Dai (2008). Magnetized self-similar 
solutions of Zhang & Dai (2008) are written as 
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where the coefficients /3r, Ptj,, Pz measure the ratio of the magnetic pressure in three directions to the gas 
pressure, i.e. Pr,4),z = Pmag,r,4>,z/Pgas- The coefficients ci, C 2 and C 3 are obtained using a set of algebraic 
equations (Zhang & Dai 2008): 
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Here, 7 and s are adiabatic index of the gas and mass loss parameter, respectively. Also, / measures the 
degree to which the flow is advection dominated. Now, we can substitute the above magnetized self-similar 
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3 ANALYSIS 

We now study the root mean radial velocity square of the clumps as a function of the radial distance 

for different values of the input parameters using the main equation (3). In all Figures, we assume coefficient 
of the viscosity is a = 0.1 and the adiabatic index is 7 = 1.4. Moreover, the mass of the central object is 
fixed at one solar mass. The coefficients of the drag force are = 5 X 10 ^ and F^ = 2.8 x 10 

Figure 1 shows the root mean radial velocity square of the clumps as a function of the radial distance 
normalized by Rsch for different values of Pr whereas the other magnetic parameters are fixed as = 
P(j, = 1. This Figure indicates that the value of increases with Pr, although its variation is not very 

significant. 

In Figure 2, we assume that the radial component of the magnetic field does not exist and the toroidal 
component is fixed, i.e. Pr = 0 and P^ = 1. We can then vary the parameter P^ to study its effect on the 
radial dynamics of the clumps. Again, we see that the value of increases as the vertical component 

of the magnetic field becomes stronger, though its variation with P^ is less significant at large values of P^. 
Moreover, at the inner parts of the system clumps are radially moving faster as the parameter Pz increases. 

Dependence of the root mean radial velocity square of the clumps on the variations of the the toroidal 
component of the magnetic field is more complicated as it has already been explored by (Khajenabi et al. 
2014) for a purely toroidal configuration. In Figure 3, we assume that /3j. = 0 and Pz = 1, but different 
values of are considered. In comparison to the previous study (Khajenabi et al. 2014), here, the vertical 
component of the magnetic field is considered too. The value of the root mean radial velocity square of the 
clumps decreases at the inner parts of the system with P^ whereas the value of increases at the outer 

parts of the system. 

Since the radial and the rotational velocities of the gas component strongly depend on the amount of 
the advected energy, then obviously the clumps experience different values of the drag force depending on 
the advection parameter /. We explore dependence of on the parameter / for different magnetic 

field configurations in Figure 4. For purely radial or toroidal magnetic field geometries, the value of 
strongly increases with the amount of the advected energy. But for purely vertical magnetic field, this trend 
changes to a reduction of the mean velocity square of the clumps with increasing the parameter Pz- However, 
this reduction is not very significant. Thus, one may conclude that the value of generally increases 

as the flow becomes more advective even in the presence of all three components of the magnetic held 
(Figure 5). 

In all of the figures, most of the curves cannot extend to the very inner parts. It is actually because of 
the limitation of similarity solutions. We described the gaseous component using self-similar solutions and 
as we know similarity solutions are valid only at the regions far from the boundaries. In other words, these 
similarity solutions for gas component are not valid at the very inner parts. But at the intermediate regions, 
similarity solutions represent dynamics of the gas flow with a very good accuracy. Since dynamics of clumps 
in our model is determined mainly due to interaction of clumps with the gas and similarity solutions for the 
gas are not valid at the inner boundary, we not investigate fate of the clumps at the inner parts based on our 
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Fig.l Root of the averaged radial velocity square {v^) of clumps versus the radial distance for a 
central object with one solar mass. Different values for the parameter /3r are considered and each 
curve is labeled by the corresponding value of this parameter. The rest of the input parameters 
are 7 = 1.4, s = —0.5, / = 0.9, /3z = 1 and /3^ = 1. 


4 CONCLUSIONS 

We studied dynamics of an ensemble of cold clumps embedded in a hot magnetized accretion flow. Although 
magnetic field has a vital role in stability and confinement of cold clouds, their role on the orbital motion 
of these clumps has not been studied. In our work, properties of the gas component is modified in the 
presence of a global magnetic field and so, the drag force on each clump changes accordingly. Comparing 
to the previous study by Khajenabi et al. (2014) who assumed the toroidal component of the magnetic 
field is dominant, we showed that both the radial and the vertical components of the magnetic field also 
lead to some changes in the averaged radial velocity square of clumps. The value of increases 

with increasing the strength of the radial and the vertical components of the magnetic field. Moreover, 
velocity dispersion of clumps increases as the flow becomes more advective when all components of the 
magnetic field are considered. Although results of our analysis are not directly applicable to the real systems 
because of limitations of this simplified model, the present study clearly demonstrates the importance of the 
magnetic held in the dynamics of clumps which can not be neglected. The results of the paper are obtained 
within the conditions of strong coupling and simplification of the magnetic field. It is also possible to relax 
these simplifying assumptions but then it would be very unlikely to obtain analytical solutions which is our 
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Fig. 2 Same as Figure 1, but for different values of j3z ■ Here, we have /3r = 0 and = 1. 

As the clumps move toward the central black hole, they will gradually accumulate at the inner parts 
because of the tidal disruption of the black hole’s gravitational field. In fact, tidal disruption determines the 
inner edge of the clumpy disc. In the presence of a global magnetic field, we find that on the average the 
clumps are radially moving faster in comparison to a similar configuration but without magnetic field. WCL 
calculated the capture rate of clumps and found that it is directly proportional to the ratio of 
Thus, presence of a global magnetic field increases the capture rate of clumps, but the level of enhancement 
depends on the detailed input parameters as we explored in this study. In other words, capturing clumps is 
faster when magnetic fields are considered. 
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Fig. 3 Same as Figure 1, but for different values of j3cj >. Here, we have fSr = 0 and /3z = 1. 
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Fig. 4 Root of the averaged radial velocity square of clumps versus the radial distance 

for a central object with one solar mass Here, we explore dependence of the root of the aver¬ 
aged radial velocity square on the amount of the advected energy for different magnetic field 
configurations. 
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Fig. 5 Same as Figure 4, but all three components of the magnetic field are considered. 
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